Cubatic phase for tetrapods 
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We investigate the phase behavior of tetrapods, hard non-convex bodies formed by 4 rods con- 
nected under tetrahedral angles. We predict that, depending on the relative lengths of the rods these 
particles can form a uniaxial nematic phase, and more surprisingly they can exhibit a cubatic phase, 
a special case of the biaxial nematic phase. These predictions may be experimentally testable, as 
experimental realizations of tetrapods have recently become available. 
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I. INTRODUCTION 

The simplest liquid crystalline phase is the nematic. It 
is a spatially homogeneous phase in which the orienta- 
tions of the non-spherical component particles, e.g. rod- 
like or disk-like colloids, are distributed in an anisotropic 
fashion. More precisely, they are oriented around a pre- 
ferred axis yielding a phase with macroscopic uniaxial 
optical anisotropy. However, this represents only the 
simplest form of rotational symmetry breaking. When 
in addition the cylindrical symmetry around the nematic 
director is broken, the phase that results is the so-called 
biaxial nematic phase. As the name biaxial suggests, 
there are now two preferred axes, which are mutually 
perpendicular. Biaxial phases can be expected if the con- 
stituent particles themselves do not have (effective) cylin- 
drical symmetry, but are only invariant under a limited 
number of discrete rotations 0, Q • Biaxial phase may 
also appear in mixtures of rod-like and disk-like particles. 
Each of the two components individually will form uni- 
axial nematic phases at sufficiently high densities. When 
mixed they will do the same, but their mutual interaction 
is such that the preferred orientation axis for the rods is 
perpendicular to that of the disks 0, 0, IU, El • 

This still does not exhaust all possibilities for spatially 
homogenous liquid-crystalline phases. Frenkel pj pro- 
posed that particles consisting of three identical rods, 
connected at right angles at their center, should form a 
stable high-density phase with cubic orientational sym- 
metry. This liquid crystalline phase is referred to as a 
cubatic phase. It is a special case of the general biaxial 
phase, since there are now three mutually perpendicular 
axes of symmetry that are equivalent. 

This model has subsequently been generalized to cross- 
like particles (" Onsager crosses" ) , in which the three rods 
can have unequal lengths |J . These particles show a sur- 
prisingly rich phase behavior. Not only do they form a 
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cubatic phase in the case that the three rods have ap- 
proximately equal lengths, but they can also show rod- 
like and plate-like behavior by forming uniaxial nematic 
phases if one, respectively two, rods are dominant in de- 
termining the shape of the particle. Moreover, at higher 
densities, these uniaxial nematic phases become unstable 
and different types of biaxial nematic phases are formed. 

Unfortunately these lower-symmetry liquid-crystalline 
phases have, thus far, not been observed in experiment. 
The main problem seems to be that, unlike the rod-like 
colloidal particles that form nematics, cross- like particles 
that are both rigid and sufficiently monodisperse could 
not be made in a sufficient quantities to allow a system- 
atic study of their phase behavior. In particular the cu- 
batic phase has not yet been observed in experiments, 
although simulations have suggested that a phase with 
this symmetry may exist in a system of disk-like parti- 
cles Q. 

Recently, however, Alivisatos et al. have reported 
the synthesis of colloidal CdTe tetrapods [ljj. These 
particles could be made with a high yield and with 
well-controlled nanoscale dimensions. The experimental 
tetrapods consists of a small crystalline body from which 
four arms grow under tetrahedral angles. Since these 
arms are also crystalline, the tetrapods are fairly rigid 
and, with a suitably chosen solvent (and proper steric 
stabilization) , should behave as rigid hard-core particles. 

In this paper we consider the liquid crystalline behav- 
ior of tetrapods. For simplicity we work in the Onsager 
limit of large aspect ratios and only take into account a 
hard-core interaction. We assume that the particles are 
monodisperse, but we treat all possible combinations of 
relative lengths for the arms of the tetrapod. We focus 
here on a bifurcation analysis, which gives us an upper 
limit to the stability of the isotropic phase and yields an 
indication of the nature of the more stable liquid crys- 
talline phases. Using this analysis, we argue that the 
tetrapods of ref. [To| should, under certain conditions, 
form cubatic phases. 

In section [H] we justify the main assumptions of the 
model and derive an expression for the Hclmholtz free 
energy. We make use of rotation matrix elements, of 
which the main properties and conventions are briefly 
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described in the appendix. In section Hill we perform the 
stability analysis of the isotropic phase and interpret the 
results, and conclude with a discussion of the main results 
in section llVl 



II. THE MODEL 

To analyze the phase behavior of hard tetrapods, we 
need expression for the free energy of this system. In 
general, this is an intractable problem. However, for 
tetrapods with sufficiently slender arms we can make the 
same assumptions that were introduced by Onsager in 
the context of the isotropic-nematic transition of thin 
hard rods [Tj. Onsager showed that for a fluid of par- 
ticles with large (strictly speaking, infinite) length-to- 
width ratio, the excess free energy can be truncated at 
the second virial coefficient level. In the case of hard-core 
interactions this is equivalent to assuming that if one ran- 
domly places particles with a given density in space, the 
probability that three particles mutually overlap is neg- 
ligibly small. 

Although we have assumed large aspect ratios for the 
arms that constitute a tetrapod, it is not immediately ob- 
vious that the second virial approximation is valid. How- 
ever, since tetrapods are essentially objects with an open 
structure and consist of four connected rod-like particles, 
one would expect that if two particles overlap with each 
other this is mainly due to a single arm of one particle 
that overlaps with a single arm of the other particle. 

The validity of this plausibility argument is confirmed 
explicitely for Onsager crosses with three equally long 
arms in the isotropic phase 12]. In a detailed analysis 
it is shown that for a length-to-width ratio of about one 
thousand the probability that, under the constraint that 
two particles overlap, more than a single pair of the arms 
are overlapping, is less then a percent and decreases for 
increasing aspect ratios. In other words, in the limit of 
large aspect ratios one can describe the particle-particle 
interaction in terms of independent pairs of arm-arm in- 
teractions only. 

What Onsager showed for elongated particles in the 
isotropic phase is that the asymptotic limit of the third 
virial coefficient B3 can be expressed in terms of the 
second virial coefficient and the aspect ratio by B3 = 
B%0(D/L\og(L/D)) [l3. This has been confirmed by 
the calculation of virial coefficients for long spherocylin- 
ders The fourth and fifth virial showed a similar 

dependence B n = B^ 1] 0(D/L). Since the interaction 
between Onsager crosses in leading order is determined 
by single rod-rod interactions the same behavior should 
be observed for Onsager crosses, as indeed is found, pro- 
vided aspect ratios are of the order thousand In 
addition, the main contribution to the higher virial co- 
efficients stems from the so-called ring diagrams, which 
would lead to same scaling behavior for higher order virial 
coefficients. Hence corrections to the free energy due to 
the simultaneous interaction of three or more particles 



is an order D/L smaller than the second virial contribu- 
tion and can therefore be neglected in the limit of infinite 
aspect ratios. 

In summary then, the assumptions that virial expan- 
sion of the free energy can be truncated at the second- 
virial level and that the interaction between particles can 
be considered as a sum of pair interactions between the 
"arms" of the particles, become exact in the limit of infi- 
nite aspect ratios. For large but finite aspect ratios, these 
assumptions should constitute excellent approximations 

The truncation of the virial series leads to the following 
free-energy functional for homogeneous systems 

/?/M=Af+ /d^(0)lnV(^) 

(1) 

+ ip J dOj J dft 2 V(fiiM^2)£(^i,0 2 ) 

Here / is the Helmholtz free energy per particle, which 
is a functional of ip the orientational distribution func- 
tion (ODF). This ODF is a measure for the fraction of 
particles with an orientation J7, which is shorthand for 
the three Euler angles (a, /? 7 ) required to specify an ar- 
bitrary orientation in a fixed reference frame, and is nor- 
malized to unity. (3 = {ksT)~ l is the inverse tempera- 
ture, p the number density and / the ideal gas term that 
does not explicitly depend on the ODF. The second term 
corresponds to the orientational entropy, while the third 
term is associated to the translational entropy through 
the kernel £(Sli, fi 2 ) describing the over space integrated 
interaction of two tetrapods with orientations f2i and f2 2 - 
The interaction of hard-core objects is taken into ac- 
count via the excluded volume £(f2i,f2 2 ). This volume, 
is defined as the volume around particle 1 (with orien- 
tation orientations Oi ) that is inaccessible to particle 2 
(with orientation orientations f2 2 ). For two slender par- 
ticles, with lengths L\ and Li and diameters D\ and 
D2 respectively the excluded volume is, to leading order, 
given by 

L 1 L 2 (D 1 + D 2 )\wij\, (2) 

where 7 is the angle between the long axes of the par- 
ticles. Corrections to this expression are of order D/L. 
Since we are mainly interested in the limit of large aspect 
ratios, we restrict ourselves to the leading order only. 

In the case of tetrapods the excluded volume is, of 
course, more complicated. First of all we now have four 
rod- like arms. Wc will assume that the arms can be ap- 
proximated by cylinders with identical diameter D, but 
possibly different lengths Li, where i = 1,2,3,4. With 
the assumption that the arms of tetrapods, whilst con- 
nected, interact independently, the excluded volume of 
two tetrapods becomes a sum over the pairwise excluded 
volumes of the arms 

^ 1 ,0 2 ) = ^2i| 1) Lf ) S|sin 7 ^|, (3) 

where the superscript refers to the particle and the sub- 
script to the arm of the tetrapod, hence 7 2 12 is the angle 
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between the ith and jth arm of particle 1 and 2 respec- 
tively. For the isotropic phase this leads to a simple ex- 
pression for the second virial coefficient B 2 - It is equal to 
half the excluded volume averaged over all orientations, 
as it is simply equal to the sum of second virials for all 
pairs of rods 

B * = ^(l>) ■ ( 4 ) 

We note that we can expand any ODF depending on 
il in terms of a linear combination of Wigner rotation 
matrix elements T> l m n (f2). In what follows, we adopt the 
convention used by Brink and Satchler in the description 
of the Wigner matrices [lj] . 

In order to make use of the free energy functional Q, 
we need to rewrite the interaction @ in terms of these 
functions. To this end, we introduce an arbitrary refer- 
ence orientation of a tetrapod. We denote the directions 
of the arms with length Li of the reference tetrapod by 
the unit- vectors ej. For simplicity we assume that e\ — z 
is the positive z-direction and e 2 is lies in the xz-plane 
with a positive a;-component. Since by definition the mu- 
tual directions are under tetrahedral angles this fixes all 
directions. 

This allows us to interpret the orientation SI of a par- 
ticle, as the one we would obtain if we take the reference 
particle and rotate it over the Euler angles denoted by 
fl. Additionally we can also interpret f2 as the actual 
rotation matrix, so the directions of the arms of a par- 
ticle become fie^. Finally, we introduce the rotations gi, 
such that &i = giZ. Note that these rotations gi are not 
uniquely defined, since only two of the Euler angles are 
required in order to satisfy the restriction. However, this 



has no effect on the final result. 

It is obvious that the excluded volume Q of two 
tctrapods cannot depend on both orientations indepen- 
dently, but only on the relative orientation fl^ 1 Q 2 . Hence 
we can rewrite the excluded volume as 

£(n7 x Q 2 ) = £(ni,fla) 
= ^L^L^D^ xn 2 ef | 

i,j (5) 
= ]>>f^|z« x g-^n 2 gfz^\. 

i,j 

With the aid of this form we can now expand 
the excluded volume in terms of rotation matrix ele- 
ments V l m n (Q), by introducing the expansion coefficients 

£(a) ee E^vljn), (6) 

where I — 0, 1, • • • , oo, and — / < ?7i,rt < I. Strictly 
speaking, the rotation matrix elements also defined for 
half-integer "spin" values. However, for reasons of sym- 
metry these can be omitted fl4|. 

The expansion coefficients can be evaluated by using 
the orthogonality relation i|A2|) for the rotation matrix 
elements 

97 + 1 f 

E hm>n = J dnv l m * n (n)s(n). (?) 

Substituting the expression © and changing the integra- 
tion variables we get 



E, 



21 + 1 
8tt 2 



2DY J L l L 3 [ dnV^g&gj 1 )^ x ilz^\ 
^ 2D J^ (E^^oJ \ T,L 3 Vl q {g 3 )\ J dfi^yfi)|sin/3| 



21 + 1 
~8 



(8) 



where we made use of the symmetry relation (JAlf) and 
closure relation l|A3|l and replaced the cross product by 
its representation in Euler angles |sin/3|. 



ducing the following shorthand notations 

E hm ee J2LiV l mfi ( gi ) (9) 

i 

2i + \ r 



dpPi(cos[3)sm 2 (3 (10) 



In order for the integral to be non-zero, it is required 
that both indices of the rotation matrix element are zero. 
This is a special case for which the function reduces to 
a Legendre polynomial V l (S1) = P;(cos(/3)). By intro- 



the expansion coefficients of the excluded volume can be 
written in a compact form as 



El, m ,n — (^E>)fiiE* m Ei :n . 



(11) 
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The integral that remains can be readily evaluated (See 
[TH Eq. 7.132.1) and is only non-zero for even value of I 

tt(4Z + 1) /2A 2 
m ~ ^ (Z + 1)(2Z - l)2 4i + 2 \l ) (12) 

m+i = o. 

Finally we introduce another shorthand notation by 
using the kernel © as a functional acting on an arbitrary 
function ip(Q) 

£[if)\(Cl) = J dQ l £(Q'~ 1 n)ip(Q / ). (13) 

In particular wc allow it to operate on a rotation matrix 
element. Using the expansion © and the properties (|A1|) 
and l|A3J) this can be manipulated to yield 

p 

Note that this generates a linear combination of rotation 
matrix elements with the same value for / and rn. In other 
words each set V l m n with n = —I, ■ ■ ■ ,1 forms a subset of 
functions that is invariant under the functional operator 
of the excluded volume. We can go one step further by 
evaluating eigenfunctions of the excluded volume. Using 
the special form of the coefficients <|llfl . one can easily 
check that for each combination of I and m at most a 
single eigenvector xin exist with a non-zero eigenvalue A; 

X l m (Sl) =J2Ei, P V l m>p (Q) (15) 
p 

h = {2D)mY,Kv E i,p- ( 16 ) 
p 

Note that the eigenvalue is independent of m and in spe- 
cial cases also might become zero as for instance for odd 
values of I. 

III. BIFURCATION ANALYSIS 

The thermodynamically stable phase of our model, is 
described by the ODF that minimizes the free energy Q). 
Usually this free energy is not known exactly and one 
uses a truncated expansion as an approximate function. 
There is however one exception: the isotropic phase. In 
the limit of infinite dilution, particles do not interact and 
therefore each orientation has the same probability, hence 
the ODF is merely a constant. 

With the aid of a bifurcation analysis we can determine 
an upper limit to the stability of the isotropic phase. 
To this end, we make an expansion of the ODF around 
the stable isotropic solution. Rather than inserting this 
into the free energy QJ, we use this to find solutions of 



the stability equation that is obtained as the functional 
derivative of the free energy with respect to the ODF 

sW){ m] - x I dQm } = ' (17) 

where A is a Lagrange multiplier to take care of the proper 
normalization of the ODF. Evaluating this expression 
and using the definition (|13(l this gives us 

\n(ip(Q)) + p£[ip](Q) = A. (18) 

For both the ODF ip and number density p we take 
the formal expansion in a small parameter e 

ij) = -0o + etpi + e 2 tp 2 + ■■■ (19) 
p = Po + epi + ( 2 p 2 + ■■■ ■ (20) 

Here we use i/'o = g^r a s the ODF for the isotropic phase. 
By inserting these expansions in the stability equation 
(|18fl and grouping terms for each power in e, we obtain 
the bifurcation equations. By solving these we can find 
the lowest density po at which a symmetry breaking mode 
exists that lead to a lower free energy than that of the 
isotropic phase. In the case of a non-zero value for p\ this 
is sign of a first order phase transition at a lower density 
and hence po is the upper limit for the meta-stability of 
the isotropic phase. 

The zeroth-order bifurcation equation merely states 
that the isotropic solution ip = tpo is a solution of the 
stationarity equation. The first order bifurcation has the 
form of an eigenvalue problem 

^ + Po £[ Vl ] = o, (21) 

Wo 

where we have already eliminated the constant contribu- 
tions. Since we are interested in a non-trivial solution 
that leads to the lowest possible positive value of po, we 
only need to consider linear combinations of eigenfunc- 
tions (|15f) . In particular we need to find the one that has 
the largest absolute value among all negative eigenvalues. 

One can show that, for the case of tetrapods, there are 
only two eigenvalues that can fulfill that requirement, 



namely the ones corresponding to I = 2 and I = 4 

A 2 = - — DL 2 {AR 2 - 1) (22) 
6 

A 4 = - I |^D£ 2 (8riR 2 + i). (23) 

For practical purposes we used here L = . Li , and 



R 2 = {Y.i L2 i)/ L2 - Since A 2 = A 4 for R 2 = ^_ an d 
by construction \ < R 2 < 1, we need to distinguish two 
types of tetrapods. The ones with R 2 > f° r which 
we need to consider modes related to I = 2 and the ones 
with R 2 < for which the modes with I = 4 are the 
important ones. For the special case R 2 = we would 
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FIG. 1: The density at which the isotropic phase becomes 
unstable with respect to different modes characterized by I 
for a subset of particles with equal volume and arm lengths 
L2 = Lz — La = (L — Li)/3. The solid line (I — 2) denotes 
the isotropic-nematic instability and the dotted line (I = 4) 
indicates the isotropic-cubatic instability. 



actually need to consider combinations of both, which 
makes the analysis somewhat more involved, but since 
this is not going to lead to new insights we will not treat 
it separately. 

As an illustration we show in Fig.^the location of the 
four lowest instabilities of the isotropic phase for a spe- 
cific class of particles, i.e. those for which L2 = L% = 
L4 = (L — Li)/3. This set includes the fully symmet- 
ric tetrapod {Li/L = 1/4) and a limiting tetrapod with 
only a single arm (Li/L — 1). Note that the volume of 
each particle is the same v — (tt/4)D 2 L, but the density 
at which the isotropic phase becomes unstable for the 
symmetric tetrapod is approximately thirty times higher 
than that of a tetrapod with a single arm. The volume 
fraction <fi — pv, however, is proportional to D/L, which 
means that in the limit of large aspect ratios the transi- 
tion takes place at small volume fractions. 

Although we now know the upper limit to the stability 
of the isotropic phase, we do not yet know which are 
the symmetry breaking modes. To find these, we need 
to perform a second order bifurcation analysis, i.e. solve 
the equation 



■po£[ih]+Pi£[ipi 



dVL-^ (24) 



employing the general solution of the first order bifurca- 
tion equation given by jUJ 



Vh(ft) 



/ / WnAm 



(25) 



where the c m are some complex constants and I is ei- 
ther 2 or 4. Substitution in the second-order bifurcation 
equation 124|) . multiplying with x* an d integrating over 



the orientation fl, gives us a set of coupled non-linear 
equations in the coefficients c m and constant p\ 



d^X l :x l n =2^2/ Mtftf. 



(26) 



We can also extract the value of pi , if we use the complete 
function in stead of Xn 



Pi 



1 / dfltpf 
2A^ J (ffitf 



(27) 



It is important to realize that there is a restriction on p\ . 
The reason is that a non-zero value of p\ is associated 
with a first order phase transition and that in order to 
follow the solution towards lower densities we need a non- 
positive value, hence p\ < 0. 

Let us now consider the case R 2 > with I = 2 
in the preceding equations. The set of equations i|26*|) 
can be solved analytically and yields only a single non- 
trivial solution, which is degenerate since all rotations of 
a solution are also solutions of the set of equations. For 
a conveniently chosen reference frame the solution is 



Vi(n) = ct>xg(n) 

It bifurcates at a reduced density 
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P0-B2 



4R 2 - 1 ' 



(28) 



(29) 



Note that in the limit of a single arm (R 2 = 1) this 
reduces to the correct result for uniaxial rod-like parti- 
cles. This particular solution is invariant under rotations 
about the z-axis, and hence one can expect that the sys- 
tem shows a phase transition from an isotropic to an 
uniaxial nematic phase. Inserting the solution in the ex- 
pression l(2*7|l for pi we obtain 



Pi 



-7T 3 DC 



Li +L 2 - L3- Li]: 



14« 

[L\ — L2 + L3 — L^][Li — L2 — L3 



(30) 



Li}. 



From this result it follows that if the sum of the lengths 
of two arms equals the sum of the length of the two re- 
maining arms we find p± = 0. These special particle 
configurations could therefore possibly lead to a continu- 
ous phase transition and be the source of a biaxial phase. 
Whether this is scenario really applies, cannot be deter- 
mined from this analysis. One could resolve this issue 
either by solving higher order bifurcation equations or 
by a full numerical minimization of the free energy. This 
falls outside the scope of the present paper. For other 
particle configurations the requirement of non-positive 
values for p\ will fix the sign of the coefficient cq, which 
can be positive or negative, and hence fully determine the 
solution (|29|l . because the magnitude will only depend on 
the choice of normalization. 

Similar to the case for Onsager crosses we can in- 
terpret this phenomenon in terms of rod-like and disk-like 



FIG. 2: Phase diagram of the nature of symmetry breaking 
modes leading to the instability of the isotropic phase. Each 
point inside the tetrahedron corresponds to a given shape of 
the tetrapod. A vertex corresponds to a single arm, the oppo- 
site plane to a tetrapod with three arms. For particles inside 
the sphere a cubatic phase is expected, outside the sphere a 
nematic phase which is either rod-like, if it contains a vertex, 
or platelet-like if it does not. Both species are separated by 
planes denoting the particles that might have a continuous 
phase transition and show biaxial behavior. 



behavior. For each of the four arms of the tetrapod, we 
can determine the nematic order parameter, which is de- 
fined as the average value of | cos 2 (9) — i with 9 the angle 
between the nematic axis and the direction of the arm. It 
can be shown that, up to a positive normalization factor, 
this is proportional to co(4Li — L)/3 for all arms. Making 
use of permutations of arms it follows that for positive 
values of cq the longest arm has the largest nematic or- 
der, while for negative cq it would be the shortest arm. A 
special limit of the former, is the case where L\ is much 
larger than the other three. It is obvious that for such 
particles this longest arm will dominate the behavior and 
the tetrapod behaves as a single rod-like particle. The 
other extreme occurs when L\ is much smaller than the 
other three. In that case the competition among those 
three arms does not allow any of them to dominate and 
there is a preference for them to be perpendicular to the 
nematic axis and, as in the case of disks, it is the smallest 
dimension that determines the orientation of the particle. 

The set of equations l|2rj|) can also be solved for the 
case of R 2 < j^, but then we need to put I — 4 in the 
equations i|25|) - (|27ll . The bifurcation density can easily 
be determined from the proper eigenvalue (|23|1 



PoB 2 



2592 



80i? 2 + 1 ' 



(31) 



But instead of having a single family of solutions that 
solve the equations (|26|) . we now have two. The first 
family corresponds again to a uniaxial nematic phase for 
which 



FIG. 3: Cross-sections of the phase diagram for constant 
value of k, from left to right, top to bottom for decreasing 
values. The dark circular areas correspond to the cubatic 
solution, the light and white areas to the two types of nematic 
solutions, being rod- and disk-like respectively. 



is the particular solution invariant under rotations about 
the z-axis, and cq < in order for p\ to be negative. 
Contrary to the previous case Co does not change sign. 

A particular solution of the second family of solutions 
is given by 



^ (O) = c \ x 4 (n) + - (xtm + xUn)) 



(33) 



(32) 



This solution is only invariant under discrete rotations 
over 7r/2 about the x-, y-, and z-axes. It therefore cor- 
responds to the cubic symmetry group and hence to a 
cubatic phase. Also in this case, the constraint on p\ 
results in a negative value for cq. 

This analysis allows us to sketch a tentative phase di- 
agram (see Figure [Sj. In this a schematic figure, we in- 
dicate the nature of the symmetry-breaking modes that 
lead to the instability of the isotropic phase. We can 
characterize the particle shape by the normalized val- 
ues U — Li/ L with the constraint l\ + li + I3 + Z4 = 1. 
Each possible particle shape, i.e. combination of rela- 
tive lengths of the arms, corresponds to a point inside a 
tetrahedron, which is the projection of the 4-dimcnsional 
"shape" space. A vertex of the tetrahedron represents 
the limit of a tetrapod with only a single arm, the tri- 
angular plane opposite to the vertex contains all particle 
shapes for which that same arm has zero length while the 
remaining three arms have non-zero lengths. In Figure 
13 we have made cross-sections corresponding to planes 
with constant value for one of the Zj. 

Nine distinct regions can be identified: a spherical re- 
gion corresponding to R 2 < ™, where there are two 
modes that lead to the instability of the isotropic phase, 
one of which yields the cubatic solution. In Figure[3]these 
regions are indicated as the dark circular areas. Outside 
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these regions, we find the particles for which there is only 
a mode with uniaxial symmetry that makes the isotropic 
phase unstable. There are four equivalent areas close 
to a vertex, shown as white in Figure |3 that represent 
tetrapods where the solution (|28|l has a positive sign, and 
four equivalent areas, indicated by light gray in Figure^! 
where the solution has a negative sign. They are sepa- 
rated by the planes(lines) leading to "biaxial" particles. 



IV. DISCUSSION 

The present work suggests that the availability of 
nanocrystalline tetrapods with a well controlled size and 
shape , may make it possible to observe cubatic liquid 
crystalline phases in experiments. Of course, the present 
" Onsager-style" analysis only becomes exact in the limit 
of very slender, rigid arms. Within this approximation, 
we have determined the upper limit to the stability of the 
isotropic phase by means of a bifurcation analysis. In ad- 
dition, we have determined the nature of the fluctuations 
that cause the instability. 

Roughly speaking we can distinguish two types of 
tetrapods, the ones for which the arms have approxi- 
mately the same lengths and the ones where one or more 
arms are significantly longer than the others. For the first 
group the isotropic phase becomes unstable with respect 
to a distortion with either nematic or cubatic symme- 
try, while for the second group only a nematic symmetry 
comes into play. 

In general the transitions will be first order since we 
obtained a non-zero value of the first order shift in density 
pi along the bifurcating solutions, indicating the presence 
of a v.d. Waals loop. An exception might be formed by 
the particles located in the planes in Figure |5J Although 
the results of the bifurcation analysis cannot guarantee 
that the symmetry of the fluctuations that lead to the 
instability will also be the symmetry of the more sta- 
ble phase, the experience in a similar study of Onsager 
crosses has shown that the bifurcation analysis has a high 
predictive value Q . In addition, it also strongly indicates 
that a system of tetrapods with approximately identical 
arms will have an isotropic to cubatic phase transition, 
which in principle could be verified by a full minimization 
of the Helmholtz free energy functional. 

Based on the results of the Onsager crosses, we pre- 
dict that at higher densities a system of tetrapods will 
undergo additional transitions to phases with yet lower 
symmetry, ultimately arriving at the phase where only 
arms with identical lengths are aligned. In particular 
three intermediate phases might appear that are invari- 
ant under 2-fold, 4-fold, and 6-fold rotations from the 
cubic group. In some cases these phases could actually 
preempt the isotropic-to-nematic or isotropic-to-cubatic 
phase transition. This is most likely to happen for par- 
ticle shapes close to the planes and/or surface that sep- 
arate the different regions in the phase diagram. A full 
numerical free-energy minimization would be required to 



confirm the existence of these phases and to ascertain 
whether a transition from isotropic phase to any of these 
four phases is possible. 

The formation of a cubatic phase even for fully sym- 
metric tetrapods may seem surprising. Naively, one 
might expect to observe a "tetrahedratic" liquid crys- 
talline phase. The reason why the latter phase does 
not appear here is related to the fact that in the On- 
sager approximation presented here the arms of a parti- 
cle effectively interact independently with those of other 
particles. Essentially we are therefore insensitive to the 
details on how the arms are connected. Although the 
relative orientation of the arms within each particle is 
fully accounted for, we cannot specify that the arms are 
connected end-on. The same results would therefore be 
obtained for particles in which the four rods were con- 
nected at their mid-points. Such particles would in fact 
have cubic symmetry and are hence are unable to form a 
phase with tetrahedral symmetry. 

The results presented here are only valid in the limit of 
infinite aspect ratios. For finite aspect ratios we expect 
two types of corrections: (i) contributions due to simulta- 
neous overlap of three or more arms and (ii) dependencies 
on the detailed construction of the particles. In practice, 
one would presumably require aspect ratios of the order 
one thousand or more in order to make these corrections 
negligible (see Ref. |8|]). Nevertheless one could expect 
that our qualitative findings remain valid even for smaller 
aspect ratios, as is the case for single rods. The fact that 
particles might not be perfectly monodisperse is proba- 
bly not a problem, since there is a rather broad range of 
particle shapes that gives rise to the cubatic instability 
of the isotropic phase. However, finite aspect ratios im- 
ply that also non-homogeneous phases should be consid- 
ered, in which the isotropic-to-cubatic transition might 
be preempted by crystallization. Finally, there may be 
kinetic limitations to the formation of cubatic phases of 
tetrapods, as the tetrapods are likely to become entan- 
gled at high densities, which could lead to kinetically 
arrested glass-like phases. 

APPENDIX A: ROTATION MATRIX ELEMENTS 

Here we list the main properties of the rotation matrix 
elements. For a more extended discussion we refer the 
reader to Q. 

A rotation = (a, (3, 7) is obtained by successive ro- 
tations of angles a, (3, and 7 about the z-, y-, and z-axis 
respectively. The invariant measure of the rotation is 
given by dfl = sm((3)dadf3d r y 1 with a, 7 £ [0, 2n] and 
/3£[Q,tt]. 

Symmetry relation 

P^(O) = P^jn- 1 ) (Al) 
Orthogonality relation 

J dnV l m ,* n ,(Q)V l ^ n (Q) = ^—^,5^,5^, (A2) 
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Closure relation 

i 

©^(Jfcfil) - E -DLpi^H-.n^) ( A3 ) 
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